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^ ■ Abstract 

O: 

We present results for VF-boson production at large transverse momentum at LHC and 
I Tevatron energies. We calculate complete next-to-leading-order (NLO) QCD corrections 

and higher-order soft-gluon corrections to the differential cross section. The soft-gluon 
(— I . contributions are resummed at next-to-next-to-leading-logarithm (NNLL) accuracy via 

0^ the two- loop soft anomalous dimensions. Both NLO and approximate next-to-next-to- 

pi I leading-order (NNLO) pT distributions are presented. Our numerical results are in good 

D . agreement with recent data from the LHC. 

^! 1 Introduction 

^ ■ The production of W bosons with large transverse momentum, pt, has been observed and ana- 
lyzed at the Tevatron over the past two decades, and significantly higher event rates have been 
observed as expected at the LHC over the past couple of years. This Standard Model process 
O I is a background to Higgs production and new physics and thus it is important to have accurate 
^ I theoretical predictions to exploit fully the large number of events at the LHC. The px distri- 
^ ■ bution falls rapidly with increasing pt, spanning several orders of magnitude over accessible 
• ^ ■ regions at hadron colliders. High-p^ W production has a clean experimental signature when 
^ ! the W decays to leptons, and solid predictions are needed to reduce uncertainties in precision 
I measurements of the W mass and decay width. The charged leptons in complementary pro- 
cesses involving Z bosons can be measured with somewhat higher resolution than the neutrino, 
but the observed event rate for W bosons at the LHC is as much as a factor of ten larger 
than that for Z bosons. Precise calculations for W production at large pr are also needed to 
identify signals of possible new physics, such as new gauge bosons, which may enhance the pt 
distribution at large pr- 

At leading order (LO) in the strong coupling as, & W boson can be produced with large 
Pt by recoiling against a single parton which decays into a jet of hadrons. The LO partonic 
processes for W production at large pt are qg — )■ Wq and qq — )■ Wg. 

The NLO corrections arise from one-loop parton processes with a virtual gluon, and real 
radiative processes with two partons in the final state. The NLO corrections to the cross section 
for W production at large pt were calculated in [H |2] where complete analytic expressions were 
provided. Numerical NLO results for production at the Tevatron were also presented in Refs. 
[1112]. The predictions are consistent with the data from the CDF [3] and DO [3] collaborations. 
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Figure 1: LO diagrams for the process qg — )■ Wq. 



The NLO corrections enhance the differential distributions in pt of the W boson and they 
reduce the factorization and renormahzation scale dependence. 

Beyond NLO, it is possible to calculate contributions from the emission of soft gluons. 
These corrections can be formally resummed and they were first calculated to next-to-leading- 
logarithm (NLL) accuracy in [5]. Approximate NNLO corrections derived from the resumma- 
tion were used in [6] and were shown to provide enhancements and a further reduction of the 
scale dependence. Numerical results were presented for the Tevatron in [B] and for the LHC 
at 14 TeV energy in Ref. [7]. In this paper we extend the resummation to next-to- next-to- 
leading- logarithm (NNLL) accuracy (see also [8]). A related study using soft-collinear effective 
theory (SCET) has recently appeared in [9]. 

In the next section we briefly review the NLO calculation and present numerical results for 
the pt distribution of the W at the LHC and the Tevatron. Section 3 discusses NNLL resum- 
mation for the soft-gluon corrections. In section 4 we derive approximate NNLO expressions 
from the NNLL resummation and we present approximate NNLO pt distributions for the W 
boson at the LHC and the Tevatron. We conclude in Section 5. 



2 NLO results 

We start with the leading-order contributions to W production at large pt with a single hard 
parton in the final state. The two contributing sub-processes are 

q{Pa) + 9{Pb) WiQ) + q{pc) 

and 

q{Pa) + q{Pb) W{Q) + g{p,). 

We define the kinematic variables s = {Pa + Pb)'^, t = {Pa — Q)"^, u = {pb — Q)'^, and S4 = 
s + t + u — Q"^. At the partonic threshold, where there is no available energy for additional 
radiation, S4 — )■ 0. The partonic cross sections are singular in this limit, but the divergences 
are integrable when averaged over the parton distributions in the colliding hadrons. 

The LO diagrams are shown in Figs. [1] and [21 The LO differential cross section for the 
qg — )■ Wq process is 
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Figure 2: LO diagrams for the process qq Wg. 



where 



qg^Wq 



s{N^ - 1) ^ Ll^//J ' (2-2) 




with the renormahzation scale, L the left-handed couplings of the W boson to the quark 
line, and / the quark flavor. Also Cp = {N^ — l)/(2A'"c) with Nc = 3 the number of colors. 
For the process qq — )■ Wg the LO result is 



where 



^qq^Wg - ^ \^hfa\ , 1^-4) 

t u tu 

The complete NLO corrections were derived in [1] [2]. The virtual corrections involve ul- 
traviolet divergences which renormalize as and make it depend on the renormahzation energy 
scale which we set to be Both the real and the virtual corrections display soft and 

coUinear divergences which arise from the masslessness of the gluons and the zero-mass approx- 
imation for the quarks. The soft divergences cancel between real and virtual processes while 
the collinear singularities are factorized in a process-independent manner and absorbed into 
factorization-scale-dependent parton distribution functions. 

The complete NLO corrections to the LO differential cross section can be written as a sum 
of two terms 

Eq ^°^73g'^^^'' = «'(/^«) ['^(^4) B{s, t, u, fin) + C{s, t, u, s,, fip)] (2.5) 

with Hp the factorization scale. The coefficient functions B and C depend on the parton flavors. 
B{s, t, u, fiji) is the sum of virtual corrections and of singular terms ~(5(s4) in the real radiative 
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W production at LHC NLO S =7 TeV 



W production at Tevatron NLO S "=1.96 TeV 
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Figure 3: Ratios of the ly-boson NLO px distribution with various choices of scale to the 
central result with scale = at the LHC at 7 TeV (left) and at the Tevatron (right). 



corrections. C{s,t,u, S4, fip) is from real emission processes away from S4 = 0. The NLO 
corrections are crucial in reducing theoretical uncertainties and thus making more meaningful 
comparisons with experimental data for W production at the Tevatron [21 S] and the LHC [10] 
at large transverse momentum. 

All numerical results presented in this paper are for the sum of and W~ differential 
cross sections. Initial-state parton densities are taken from MSTW2008 [11]. We use the NLO 
central sets with the QCD coupling evolved at NLO in this section. In section 4, where we 
include the NNLL logarithmic corrections, we employ the NNLO central sets with the QCD 
coupling evolved at NNLO. The W is on shell and final-state partons are integrated over the 
full phase space. 

We begin with results for W production at the LHC at center-of-mass energy 
and with pt in the range 20-500 GeV. At lower values of pt the fixed order NLO estimates 
become unreliable due to Sudakov logarithms and the comparison with experiment needs to 
include nonperturbative pt smearing. We set the factorization and renormalization scales equal 
to each other and denote this common scale by fi. In the left plot of Fig. |3]we plot ratios for 
the W-boson px distribution at the LHC with various choices of scale to the central result with 
scale /i = Pt- We display the scale variation of the NLO result with scale choices Pt/2 and 
2pT- We also show results for the choice of scale fx = \Jp^ + rriyy but note that the numbers for 
this choice of scale are very similar to those for fi = p^- We see that the scale variation is of 
the order of ±10%. The results for these ratios are almost identical at 8 TeV energy and very 
similar at 14 TeV. 

In the right plot of Fig. [3] we show results for W production at the Tevatron at 1.96 TeV 
energy. Again, we display the scale variation of the NLO result with scale choices Pt/2 and 

2pT and also y^pfH-rn^. We note that the published CDF and DO analyses [Sj H] involve older 
Run I data. The results presented here are also applicable to the higher luminosity and energy 
data from Run II. 
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Figure 4: One-loop eikonal diagrams for qg — )■ Wq. 
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Figure 5: One-loop eikonal diagrams for qq — )■ Wg. 



We will say more about the scale variation in Section 4 when we include the NNLO soft- 
gluon corrections. Another source of uncertainty in the differential distributions comes from 
the parton distribution functions (PDF). We will return to the topic of PDF uncertainties when 
we present the approximate NNLO results in Section 4.3. 



3 NNLL Resummation 



Near partonic threshold the corrections from soft-gluon emissions are dominant. These correc- 
tions can be resummed to all orders using renormalization group arguments. The resummed 
cross section is derived in Mellin moment space, with N the moment variable conjugate to S4, 
and is given formally by 



exp 



exp 



E'AN'] 



exp 



(3.1) 



where the first exponential resums the coUinear and soft-gluon radiation from the inital-state 
partons; the second exponential resums corresponding terms from the final state; the third 
exponential controls the factorization scale dependence of the cross section via the parton- 
density anomalous dimension; H is the hard-scattering function; and S is the soft-gluon function 
describing noncollinear soft gluon emission whose evolution is controlled by the soft anomalous 
dimension F^. The first three exponentials in Eq. (13.11) are independent of the color structure 
of the hard scattering and thus universal [121 [13], while the functions if, and F^ are process- 
specific [51 [T3|. More details of the resummation formalism have been given before (see e.g. 
Refs. [51 [m [U]) and will not be repeated here. 
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Figure 6: Two- loop eikonal diagrams for qg Wq and qq — )■ Wg involving the two incoming 
partons. There are two additional sets of 12 diagrams each with the same topologies that 
involve one incoming and one outgoing parton. 



We expand the process-specific soft anomalous dimensions Ts in the strong coupling as 



(3.2) 



The one-loop results, r^^\ are obtained from the ultraviolet poles in dimensional regularization 
of one-loop eikonal diagrams involving the colored particles in the partonic processes, Figs. H] 
and[5|, and were first derived in |5]. We determine the two-loop results, F^^"*, from the ultraviolet 
poles of two-loop dimensionally regularized integrals for eikonal diagrams shown in Fig. E] and 
related graphs involving other combinations of the eikonal lines (see also [8]). 
For qg — Wq the one-loop soft anomalous dimension is 



and the two-loop soft anomalous dimension is 



r(2) 

S,qg^Wq 



K 



2 S,qg^Wq ' 



(3.3) 



(3.4) 



where K = C^(67/18 — (2) — 5^//9 [16] with Ca = 3 and nj the number of light quark flavors. 
For qq — Wg the corresponding results are 



and 



S,qq^Wg 



r(2) 

S, qq^Wg 



2 5", qq^Wg 



(3.5) 



(3.6) 
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We note that the proportionahty of the two-loop soft anomalous dimension to the one-loop 
result is anticipated on general grounds from the work in Ref. [T7] (see also [HI [191 120]). 



4 NNLO approximate results 

By expanding the resummed cross section, Eq. (13. ip . in the strong coupling we derive approx- 
imate fixed-order results. In this section we present the analytical expressions for the NNLO 
expansion and use them to present approximate NNLO results for the VT-boson transverse 
momentum distribution at the LHC and the Tevatron. 



4.1 qg — > Wq 



We can write the NLO soft and virtual corrections for qg — )■ Wq as 



Q- 



(1) 

qg^Wq 

d^Q 



qg^Wq' 



TT 



of 



ln(g4/p|) 

S4 



1 

Si 



(4.1) 



The NLO coefficients c^^ and c^^ of the soft-gluon terms in Eq. (14.11) can be derived from the 
expansion of the resummed cross section and are given hj c^^ = Cp + 20 a and 
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(4.2) 



For later use in the NNLO results we also define the scale-independent part of c^f as T^^ 
-{?,/A)CF-CAH'tu/spl). 

The coefficient of the 5{s^ terms in Eq. (14. ip is given by 



r.19 



1 r'^ 

^ + Bfuf + Cr + Cfuf] + ^ In^ 



+ cf In 



(4.3) 



with Bi^, i?2^, C*!^, and Cf^ as given in the Appendix of the first paper in Ref. [2] but without 
the renormalization counterterms and using = hi{A/Q'^) = [note that the terms not 
multiplying A'^^ in Eq. (A4) for Bf^ of Ref. [2] should have the opposite sign than shown in 
that paper]. 

The NNLO expansion of the resummed cross section involves logarithms ln^(s4/py) with 
k = 0,1,2,3. The terms with k = 1,2,3 were already provided in Refs. [5l [6] from NLL 
resummation. Terms for k = were also given in Eq. (3.8) of Ref. [6] but they were incomplete. 
Now that we have expressions for NNLL resummation we can provide the complete result for 
the k = terms. 

The NNLO soft-gluon corrections for qg — > Wq can be written as 



E, 



da 



Q- 



(2) 

qg^Wq 



d^Q 



qg^Wq' 



■ a' 



qg^Wq 



(4.4) 
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qg-^Wq 
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-S4- 



with Pq = {IICa — 2nj)/3 and where we have used the two-loop constants (cf. [22]) 



101 11 7^ , ^ 
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= (CA/C^)Df , and 
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(4.7) 



Note that the difference from Eq. (3.8) of Ref. [6] is in the [l/s4]+ terms. The higher 
powers of the logarithms are the same. Also note that one can determine the scale-dependent 
6(34) terms at NNLO. These terms are also given in Ref. [6] and will not be repeated here. 



4.2 qq — > Wg 

For the process qq — )■ Wg we can write the NLO soft and virtual corrections as 
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Here the NLO coefficients cf and cf of the soft-gluon terms are cf = ACp — Ca and 

tu 
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(4.9) 



For later use in the NNLO results we also define the scale-independent part of cf as Tj'' 
-(2C^-CA)ln(tn/sp2.)-/3o/4. 



The coefficient of the 5{s4) terms in Eq. (14. 8p is given by 
1 



2Ai'i 



Bf + Cf + {Bf + D^^J) 



n 



<i<i 
C3 



(4.10) 



with Bf"^, B2^, Cf^, and D^J as given in the Appendix of Ref. [2] but without the renormal- 
ization counterterms and using /a = 0. 

Again, the NNLO expansion involves logarithms ^^^(54/^^) with k = 0,1,2,3. The terms 
with k = 1,2,3 were already provided in Refs. [3, E] from NLL resummation, but the terms 
given for = in Eq. (3.19) in [B] were incomplete. With NNLL resummation we can now 
provide the complete result. 

The NNLO soft-gluon corrections for qq — > Wg can be written as 
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where 



B^^ = C\ 



1025 



7C3 + 



79 



(4.11) 



(4.12) 



(4.13) 



432 4' / 108 ■' 8 108 

The scale-dependent ^(34) terms at NNLO were also provided in Eq. (3.19) of ^ and we 
will not repeat them here. 



4.3 Numerical results 

We begin with results for W production at the LHC at 7 TeV energy. In Fig. [7] we plot the 
M^-boson "px distribution, da/dpT- In the left plot, we compare the NLO and the approximate 
NNLO results at the LHC at 7 TeV energy with ^ = px- We also compare our results to recent 
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W-boson p distribution at LHC S =7 TeV ^=p 
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Figure 7: VT-boson approximate NNLO pt distribution at the LHC at 7 TeV compared with 
ATLAS data. 



data from ATLAS |T0]. It is evident that the effect of the NNLO soft-gluon corrections grows 
with Pt as one would expect, since the kinematical region near partonic threshold becomes 
more important at higher px- The inset plot shows that the ratio of the approximate NNLO to 
the full NLO (i.e. the K factor) grows with px-, and the NNLO soft-gluon corrections provide 
nearly a 60% enhancement at pr = 500 GeV. Since the ATLAS data use acceptance cuts 
and are normalized by the total fiducial cross section, we have to correct for these factors to 
extrapolate the experimental results for direct comparison to our pt distribution. We use the 
procedure described in [23]. We multiply the normalized ATLAS results by the total fiducial 
cross section and divide by the acceptances. It is clear from the comparison that the data are in 
very good agreement with our NNLO approximate result, which provides a better description 
than NLO alone. The ATLAS data only go up to a p^ of 300 GeV and it will be interesting to 
see data from the LHC at even higher p^- 

In the right plot of Fig. [7] we show ratios of the approximate NNLO result with the variation 
of the result with scale /i = pt/2 and 2pT relative to fi = pt- We also show the ratio with 
/i = \Jpx + f^w note that the results for this choice of scale are very similar to those for 
/i = Pt- Also, comparing the left plot of Fig. [3] with the right plot of Fig. [7] it is seen that 
the scale dependence at approximate NNLO is smaller than at NLO at intermediate values of 
Pti but it grows larger at very high pt values where the overall soft-gluon contribution is also 
larger. 

We note that results for 7 TeV LHC energy have also been provided in Ref. [9] at NNLL 
in the SCET formalism. It is not possible to do an analytical comparison with |9] because no 
analytical results are provided there. It is also important to note that NNLL means different 
things in different approaches and also if one uses different variables within the same approach. 
This has been clearly explained in the context of top quark production in the review paper 
of Ref. [23] and applies here as well. Therefore our NNLL Mellin-space resummation is not 
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Pj (GeV) (GeV) 



Figure 8: W^-boson approximate NNLO px distribution at the LHC at 8 TeV. 



equivalent to the NNLL SCET resummation of p]. Related studies for top quark cross sections 
show that different formalisms can give very different numerical results [24]. 

Furthermore, it is difficult to make a meaningful numerical comparison because most of the 
figures in [9] plot results versus parameters that are only defined within the SCET formalism. 
Non-graphical results are provided in Table 1 of |9], which provides absolute integrated cross 
section estimates for px > 200 GeV. The LO and NLO fixed-order results are of course in 
agreement with ours, but exponentiated logarithmic corrections will obviously diverge in this 
bin, thus making the comparison not very meaningful. At a purely numerical level, however, 
we note that our approximate NNLO integrated cross section at 7 TeV LHC energy above a 
Pt of 200 GeV is about 20% larger than that provided in Table 1 of Ref. [9]. 

Finally it is important to note that here we use an NNLO expansion of the resummed 
expression to obtain numerical results, in order to avoid prescriptions needed to invert from 
moment space to momentum space. Even within a given formalism such as SCET or Mellin 
moments, there is a numerical difference between using NNLO expansions and attempting a full 
resummation (again see [2lj and references therein for more details). However, this difference 
is typically smaller than the differences between different formalisms or prescriptions, so this is 
usually a point of relatively minor significance. For example both the values and uncertainties 
in [H] are very similar for the resummed and NNLO approximate results. 

We continue with the presentation of our numerical results for the ly-boson pj- distribution 
at the LHC. In Fig. [8] we show results for the LHC at 8 TeV energy. Although the overall 
Pt distribution is enhanced at 8 TeV relative to 7 TeV, the result for the ratio of approximate 
NNLO over NLO shown in the inset plot is very similar (almost identical) to that at 7 TeV. 
Also the scale dependence at 8 TeV is almost the same as at 7 TeV as can be seen by comparing 
the right plots of Figs. [7] and El 

In Fig. [9] we show the corresponding results for the LHC at 14 TeV energy. Again, the 
increase due to the NNLO soft-gluon corrections is more prominent at very high p^, reaching 
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Figure 9: VT-boson approximate NNLO pt distribution at the LHC at 14 TeV. 



90% enhancement at pt = 1000 GeV. The scale variation at the LHC at 14 TeV energy is 
shown on the right plot. At the very highest pt the scale uncertainty is significant. 

Finally, we present results for the Tevatron. The left plot of Fig. [10] shows the NLO and 
NNLO approximate px distributions at the Tevatron ai fi = pt, with the inset plot displaying 
the K factor. The right plot of Fig. [10] shows the results for scale variation at the Tevatron, 
where there is a reduction of scale dependence at NNLO relative to that at NLO in the pt 
range shown, as can be seen by comparing the right plots of Figs. [3] and [THJ 

In addition to the scale variation shown in the previous figures, one can also attempt an 
independent variation of the factorization and renormalization scales; however, this does not 
affect the range of the overall uncertainty significantly, if at all. 

As already mentioned in Section 2, in addition to scale dependence another important 
source of uncertainty comes from the PDF. Here we use the PDF sets and procedure provided 
by MSTW2008 [llj to calculate the PDF uncertainties for W production. In the left plot of 
Fig. [11] we compare the PDF uncertainty with the scale uncertainty at NNLO and also at NLO 
at the Tevatron. The scale ratios are for fi = pj'/2 and 2pT and are the same as already shown 
in the plots of the previous Tevatron figures, but displaying them together and with the PDF 
ratios makes the comparison of all the uncertainties easier. We note that for the Tevatron the 
PDF uncertainties are smaller than the NLO scale variation, but they are larger than the scale 
variation at approximate NNLO for most pt values. As mentioned earlier, the NNLO scale 
variation is consistently smaller than that at NLO. 

In the right plot of Fig. [11] we compare the PDF uncertainty with the scale uncertainty at 
NNLO and also at NLO at the LHC at 8 TeV energy (the results at 7 TeV are practically the 
same). Here the situation is somewhat different from the Tevatron in that the upper range of 
the PDF uncertainty is larger than both NLO and NNLO scale variation, but the lower range is 
smaller than both for most pt values. Also the scale variation at approximate NNLO is smaller 
than that at NLO for moderate pt values but becomes larger, particularly at the lower end, at 
very high pt- 
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Figure 10: ly-boson approximate NNLO pr distribution at the Tevatron. 



W production at Tevatron NNLO approx S "=1.96 TeV W production at LHC NNLO approx S "=8 TeV 
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Figure 11: PDF and scale uncertainties for the l^-boson pt distribution at the Tevatron (left) 
and at the LHC at 8 TeV (right). 
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The growing NNLO contribution relative to NLO with increasing pt at both Tevatron and 
LHC energies does not appear to be an artifact of perturbation theory. The extended pt range 
made accessible by the LHC brings into play a high multiplicity of hard jets in various inclusive 
cross sections, and this is expected to affect the convergence of perturbative predictions. The 
recent ATLAS data (TU] and its agreement with our theoretical results confirms this expectation. 
An analysis of the Z distribution is underway and will be reported elsewhere. 

5 Conclusions 

The transverse momentum distribution of the W boson receives large QCD corrections. Com- 
plete NLO calculations have been used in this paper to provide numerical results at LHC and 
Tevatron energies. In addition, NNLL resummation of soft-gluon corrections has been derived 
using two-loop soft anomalous dimensions. Approximate NNLO analytical expressions have 
been derived from the resummed cross section and employed to produce numerical results. The 
NNLO soft-gluon corrections reduce the NLO scale dependence at low and intermediate p^ 
where the bulk of the data is located. In the very high pt region at LHC energies the soft log- 
arithms seem to become sensitive to scale. It is possible that including hard NNLO corrections 
will improve this situation. The experimental bins at large pt are larger because the number 
of expected events and the px resolution both decrease dramatically. Recent ATLAS data (TU] 
are in good agreement with our numerical results. The higher-order results presented in this 
paper strengthen theoretical predictions for W production at hadron colliders and may prove 
significant for new physics searches in the tail of the pr spectrum. 
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